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Abstract 

We show how conformal invariance predicts the functional form of two-point 
correlators in one-dimensional periodic quantum systems. Numerical evidence 
for this functional form in a wide class of models — including long-ranged ones 

- is given and it is shown how this may be used to significantly speed up 

calculations of critical exponents. 
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Conformal invariance has been shown to be remarkably powerful in predicting the uni- 
versal long-distance behavior of low-dimensional field theories and statistical systems both 
at the classical and the quantum leveH Independent verification of conformal predictions 
by finite-size studies was made possible by Cardy's@ use of conformal mappings to relate a 
problem in one geometry to one in a finite-sized geometry. The numerical success of this 
approach offers a convincing validation of the principle of conformal invariance at critical 
pointJii. 

The main ideas leading to the conformal finite-size approach in two- dimensions may be 
summarized as follows. Under a conformal mapping z \— > w(z), the two-point correlation 
function of a scalar scaling operator <p(z) at a critical point will transform covariantly, i.e. 

((f)(z 1 ,Z 1 ),(j)(z2 1 Z2)) = \w'(z 1 )\ X \w'(z 2 )\ X ((t)(w 1 ,W 1 ), (j)(w 2 ,W 2 )), 

where x is the scaling dimension of 0. A mapping of the entire z-plane onto the surface of 
a cylinder of circumference L is achieved by choosing the particular conformal mappingi 

w = — ln(z) 

We can then relate the two-point function of the scalar operator <fi(z) in the infinite plane, 
to the two-point function of the operator <p(w), now evaluated in the cylinder geometry. 
Putting w = u + iv, so that u measures distance along the infinite cylinder and v is the 
periodic coordinate across, we may explicitly calculate the two-point function to be@ 

(2tt/L) 2x 

(<f>(ui, vi), 4>(u 2 , v 2 )) = -r — — — -rm (1) 



2 cosh ^f{ui — u 2 ) — 2 cos ^f{vi - v 2 
For large separation along the cylinder axis, i.e. assuming u\ — u 2 ^> L, the correlator 
exhibits an exponential decay 

/ 2,71 \ ^ x 
(0Ol,fl),0(w 2 ,f2)) = f-jT-j 

and we may deduce the correlation length 



e ~(2-Kx/L)(u 1 -u 2 ) 



For one-dimensional quantum systems, we now reinterpret the coordinate u as running 
along the infinite chain of sites and the width L in the v direction as corresponding to 
the inverse temperature f3, i.e. the above reviewed 2D results may be applied also to an 
infinitely long quantum chain at finite temperature. This then opens the possibility of 
checking the predictions for the critical exponents as obtained by (Q), by directly studying 
the long-distance behavior of the correlation functions. 

Analytic results for correlation functions, however, can only be given for a rather select 
group of mo delsfM and thus it is at this point that numerical studies and simulations 
become indispensable tools for the verification of conformal results. Various methods for 
the numerically challenging problem of calculating correlation functions have been devised 
over the years0, but due to the complexity of the problem, computations are still limited 
to lattice sizes of about N < 100. The extrapolation of the critical exponents from these 
finite-size data usually involves taking the thermodynamic limit by holding the density d 
fixed as N — > oo, as a first step. However, the correlation functions still exhibit an oscillatory 
behavior!'!. For the one-particle reduced density matrix p(r), this behavior leads to model- 
specific singularities in the momentum distribution n(k). Thus, the next step in extracting 
the long-distance behavior should be to let the average inter-site (or inter-particle) spacing 
1/d go to zero. The oscillations will then be pushed towards the origin. (For small lattice 
sizes or calculations at fixed density0Hil, where this limit is not practical, one can instead 
consider a suitable averaging procedure, leading to equivalent results.) 

There is, however, another interpretation of the two-point function ([I]). Most of the 
above numerical studies are performed not on a long chain with free boundary conditions, 
but on a ring of length N with periodic boundary conditions. Thus it seems natural, to use 
the periodic coordinate v to label the sites and the coordinate u to measure the time. In 
this picture, we just 'roll' the cylinder in a different way, i.e. we have a ring of length L = N 
and infinite extent in time. The equal-time correlator may then be constructed from ([!]) by 
choosing Am = u\ — u 2 = and defining Av = v\ — v 2 = r. The two-point function of ([!]) 
now reads as 



where we have suppressed the u±,U2 arguments. Thus we arrive at a prediction, not only of 
the critical exponents of a given theory, but also of the analytic behavior of the correlation 
functions for a finite-sized system on a ring. 

Before we present numerical evidence for the above functional form, let us examine in 
what limit this behavior is expected to occur. Because of the non-universal oscillations, we 
have true agreement only for 1/d —>■ 0. This limit can easily be achieved by plotting the 
two-point function as a function of r/L on a fixed scale from —1/2 to +1/2. We then expect 
the universal law (|]) to dominate as L — > oo. (Note that for fermionic systems, we use the 
envelope of the two-point function.) 

The consequences of (|3]) are impressive. First of all, we have a much stronger test of 
conformal invariance. After all, not only must the critical exponents agree, but also the 
complete analytic behavior. Let us now test our scaling law for the reduced one-particle 
density matrix p(r) and the spin-correlation (SoS r ) in a variety of models. For the s = 
1/2 Heisenberg antiferromagnet, it has been shown via bosonization techniques^ that the 
singular behavior of (SoS r ), i.e. 

(S S r ) oc (-l)V-" (4) 

is described by rj = 1. Recently, long standing conjectures of possible logarithmic corrections 
have been verified by including marginal operators in the above reviewed conformal finite- 
size analysis^!. In Fig. [I] we show a comparison of recent Monte-Carlo simulationsS^El and 
the scaling law (^). Note that the lattice sizes are still too small to completely push the 
oscillations to the origin. The agreement is, however, quite good. 

For hard-core bosons (the 5-function gas of Ref. |T6| ), the one-particle density matrix 
p(r) may be either deduced in terms of generalizations of the Painleve transcendent!!!, by 
application of the quantum inverse scattering method^ using the equivalence with the non- 
linear Schrodinger equation, or by random matrix theoryEl. In all methods the critical 



exponent is found to be rj = 1/2. Data taken from Ref. |9] again shows good agreement with 
the above scaling law in Fig. |^. 

Finally, the long-ranged g/r 2 models have been shown to have a set of critical expo- 
nents compatible with conformal predictions, both for their discrete versions such as the 
Haldane-Shastry mo delEM and for the continuum system!! - !!. The critical exponents vary 
continuously with interaction strength g and r\ g is independent of density d. In Fig. ||| we 
plot the results of Ref. |25| for p{r) and the scaling curve corresponding to @). Note that the 
r]2 = 1 curve corresponds to the Haldane-Shastry spin-chain. In addition, the continuum 
data coincides with the discrete data, the first being the low-density limit of the second. 

Another important consequence following from (0) is that Monte-Carlo simulations to 
determine critical exponents can now be restricted to a single calculation for each lattice 
size L. This is done most conveniently by studying the behavior half way around the chain. 
So, using as an example the spin-correlator in the s = 1/2 antiferromagnet, 

(S S L/2 ) oc C^X . (5) 



L 

Taking data from Kubo et al0, we show a plot of (S Sl/ 2 ) versus 1/L in Fig. [|. A least- 
squares fit gives r] = .94787 which is a reasonable fit for just 5 data points, especially in 
light of the logarithmic corrections. Note that the errorbars are estimates of the oscillation 
amplitudes and do not represent the statistical Monte-Carlo errors. Incidentally, this method 
has been used by the present authors to compute the critical exponents in the g/r 2 models 
to high accuracy!!. 

Let us close this short note with a speculation. Using our interpretation of equation (|J), 
we may compute the time correlators at a single site by choosing u\ — u 2 = t, v\ — v 2 = 0. 
Surpressing the V\,v 2 arguments, the result is given by 

'2tT 



(0(*),0(O)) = (y) 



2x / i \ 2x 



^ sinh ^ / 



By choosing v i — v 2 = L/2, the same reasoning gives the large distance asymptotics of the 
time correlator, 



2vr^ 2x 

T 



2.r 



In Fig. [5] we plot the complete behavior of (0) for hard-core bosons in the (t, r)-plane. It 
will be very interesting to compare these predictions to the calculations for the hard-core 
bose gas of Ref. [T7] and for the cr-model of Ref. EM 
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FIGURES 

FIG. 1. The scaled and normalized spin-correlator for N = 24, 30, 40. The line is the predicted 
scaling behavior for rj = 1. Note that even the nearest-neighbor correlations are correctly predicted. 

FIG. 2. The lattice version of the one-particle density matrix expansion p(r) for hard core 
bosons is compared to data from a direct calculation. 

FIG. 3. The scaled and normalized one-particle density function for the g/r 2 model is plotted 
with the predicted scaling curves. From left to right, the curves correspond to rj = 5/4 (g = —1/2 
fermionic, g = 4 fermionic), rj = 1 (g = fermionic, g = 4 bosonic), rj = 1/2 (g = bosonic), 
rj = 1/4 (g = —1/2 bosonic) and are symmetric about the y-axis. 

FIG. 4. We plot the spin-correlator halfway around the lattice, i.e. {SoSn/ 2 )i as a function 
of the inverse lattice length l/N. Data for N = 6,12,24,32,40 is taken from the Monte-Carlo 
calculation of Kubo et al. The solid curve is a fit for the Luther-Peschel r~ v result which gives the 
estimate rj = .94787. 

FIG. 5. The scaled correlation function of hard-core bosons (rj = 1/2) in the (t, z)-plane. The 
value att = 0,r = L/2is normalized to 1. 
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